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Abstract 

Basic notions regarding classical integrable systems are reviewed. An algebraic descrip- 
tion of the classical integrable models together with the zero curvature condition description 
is presented. The classical r-matrix approach for discrete and continuum classical integrable 
models is introduced. Using this framework the associated classical integrals of motion and 
the corresponding Lax pair are extracted based on algebraic considerations. Our attention 
is restricted to classical discrete and continuum integrable systems with periodic boundary 
conditions. Typical examples of discrete (Toda chain, discrete NLS model) and contin- 
uum integrable models (NLS, sine-Gordon models and affine Toda field theories) are also 
discussed. 



*Tliis article is based on a series of lectures presented at the University of 
Oldenburg in July 2011. 
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1 Introduction 



This is a pedagogical review aiming at introducing basic notions regarding classical integrable 
models. It is based on a series of Lectures given at the University of Oldenburg. These 
lectures are mainly addressed to graduate students or other researchers who wish to acquire 
some fundamental knowledge on certain topics within classical integrability [1]-[10]. 

The notion of complete integrability goes back to the celebrated Liouville's theorem for 
a system with a finite number of degrees of freedom [11]. In 70's the notion of integrability 
was extended to models with infinite degrees of freedom, that is 1+1 integrable field theories 
with generic non-linear interactions; such typical examples, which will be reviewed here are 
the non-hnear Schrodinger model (NLS), the sine-Gordon model and its generalizations, the 
affine Toda field theories (ATFT), among others (see e.g. [1, 3] and references therein). The 
idea of quantum integrability is relevant to the so called quantum inverse scattering method, 
an elegant algebraic technique mainly developed by the St. Petersburg group (sec e.g. 
[12]-[18], for a recent review see [19] and references therein). Perhaps the greatest appeal 
of integrability is that it offers an exact framework for the investigation of a wide range 
of physical systems without resort to perturbative methods, which are the most common 
techniques when studying any physical system. In addition to the significance of these models 
per se, and their rich mathematical structure there are various immediate applications and 
connections with other related areas of research, such as string and gauge theories, condensed 
matter physics, gravity etc. 

The main focus in this presentation lies on the description of classical discrete and con- 
tinuum integrable systems via the Hamiltonian formalism based on the classical r-matrix 
approach [20, 21]. More precisely, the outhne of this article is as follows: in the next Sec- 
tion some preliminary notions such as the Poisson structure are introduced, and Liouville's 
theorem is briefly reviewed. In Section 3 linear Poisson algebraic structures are introduced 
and integrable models ruled by hnear algebras, such as the Toda chain are described. In 
particular, the first integrals on motion are extracted and the time components of the corre- 
sponding Lax pairs are explicitly constructed from first principles. In section 4 we introduce 
the notion of quadratic Poisson structure, and through this notion we construct discrete one- 
dimensional discrete integrable models. An explicit construction of the conserved quantities 
and the associate Lax pairs is presented as well. In the next Section typical examples of 
physical systems ruled by quadratic classical algebras, such as the discrete NLS model as 
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well as the Toda chain (via the "dual" quadratic description due to Sklyanin), are presented. 
In Section 6 we describe classical continuum integrable models in the spirit presented in the 
previous sections, based again on quadratic Poisson structures. Typical examples such as 
the NLS and sine-Gordon models, affine Toda field theories (ATFT) and the Liouville model 
are reviewed in Section 7. In the last section a systematic prescription for obtaining suitable 
continuum limits of discrete integrable models in such a way that integrability is preserved is 
outlined. A considerable number of relevant exercises are provided for the interested reader 
throughout the text. 

2 Preliminaries 

The notion of Poisson bracket, which is used to define a Poisson algebra (see e.g. [1, 3, 11] 
and references therein) is briefiy reviewed. Consider a phase space 02n endowed with 
Poisson structure { , } and local canonical variables Qi, Pi. Given any two functions 
9{{Qi}: {Pi}) ^hc Canonical variables associated to the Poisson bracket 
takes by definition the form: 



For any three functions /, g, h oi the canonical variables the so called Jacobi identity 
holds, that is: 




(2.1) 



and clearly 




(2.2) 



{/, {g, h]] + [g, [h, f}} + {h, {/, g}} 



0. 



(2.3) 



Also, Leibniz's law, known as Poisson property applies: 




(2.4) 



Hamilton's equations of motion may be derived in terms of the Poisson brackets. Let a 
function f{{qi}, {Pi}, t), then 



But 



dH . dH 

Pi = 2.6 
dpi dqi 



(the dot denotes derivative with respect to time) therefore we conclude: 



For / being a constant of motion one then has 

/, i/} + ^ = 0. (2.8) 

In order for a system to be integrable all the integrals of motion should be in mutual 
involution, i.e. to Poisson commute. More precisely, let us recall the celebrated Liouville's 
theorem (sec e.g. [3, 11]): consider a system with n degrees of freedom in the phase space 
&2n with coordinates pi...,Pn, let also H{{qi}, {pi}) be the Hamiltonian of 

the system. The system is completely integrable if there exist n independent functions 
Ii{{qi}, {Pi}) defined on the phase space such that: 

=0, {/,, /,}=0, (2.9) 

li are in fact charges in involution and they are known as integrals of motions or conserved 
quantities given that they satisfy: 

dL, 
dt 

In this case it is possible to introduce a complementary set of functions (pi{{qj}, {pj}), 
i, j G {1, . . . ,n} [11], and we may then describe the change of coordinates on the phase 
space: 

{qi, P^) ^ (2-11) 

The new coordinates ^/j, (pij define the so-called action-angle variables, and the equations 
of motion become: 

^ = /jj = =^ Ij{t) = constant 

^ = {h, ^j} = UJ, ^ =uj,t + ^,(0) (2.12) 

for certain functions ujj{{Ii}). The time evolution of an integrable system is linear if we 
choose a suitable set of coordinates (angle-action). Although the existence of (jj, (pA is 



[h, /,} = 0. (2.10) 



guaranteed from Liouville's theorem (see e.g. [11]) their exphcit construction is an intriguing 
task depending on the specifics of the integrable system under consideration 

Having introduced the basic notions on the Poisson structures we are now ready to pro- 
ceed with the description of the algebraic setting underlying classical discrete and continuous 
integrable systems. 



3 Linear Poisson structure 

We start our review considering integrable models that obey linear algebraic Poisson struc- 
tures. A typical example in this category is the Toda chain. We first introduce the Lax pair 
[6] associated to such a system. It is a pair of elements of a Lie algebra Q, {L, A) satisfying 
the so called zero curvature condition: 



dL 



L, A 



(3.1) 



Assuming ^ is a loop algebra Q ® C(A), the associated spectral problem is then expressed as 

L(A) ip — u ip 

det (|l(A) - = 0. (3.2) 

The integrals of motion associated to such a system are then naturally provided by the traces 
of powers of L i.e. (see also [3, 22] and references therein), 

iW(A) = ir (L"(A)). 



(3.3) 



Expansion in powers of A (or j) is expected to provide all required integrals of motion of 
the system under consideration. 

It is quite straightforward to show that the latter provide a family of mutually Poisson 
quantities as long as L satisfies the following linear algebraic relations: 



|l„(A), L;,(/x)| = TabiX-ll), La{X) + Lbi/l) 



(3.4) 



L is in general an n x n matrix with entries being fields, representations of the underlying 
classical algebra, and r is an x matrix which satisfies the classical Yang-Baxter equation 
[21] 



ri2(Ai - A2), ri3(Ai) +r23(A2) + ri3(Ai), r23(A2) 



= 0. 



(3.5) 
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The r-matrix acts on the tensor product V ®V. In general the notation that will be used 
henceforth is: Ai = A (g) I, B2 — B and Ai B2 — A iSi B . Moreover, it is clear that the 
classical Yang-Baxter equation acts on V (8) V (g) V and ri2 = r 1 i.e. the r matrix acts non 
trivially on the first two spaces, r23 = I ® r and so on. 



Proposition 3.1. Based on the linear relations (3.4) show that: 

{t(")(A), t(™n/^)} = 0. (3.6) 
Proof. We recall the definition of i*^") and substitute: 

tra,{L:{X), L-(/.)} = nm ir„, (Lr^(A){L„(A), L,{fi)] L^-\fi)) = . . . (3.7) 
The latter expression becomes after recalling (3.4) 

... = nmtrab[L:-\X)[rab{\-fi), L,{X) + L,{fi)]L^-\li)^ 
= nm trab[L:-\X)L^-\fj,)rab{X - /x)(L„(A) + L,{fj,))) 

- nm tr,,(^L:-\X)L'^-\fj,){L,{X) + L,{fj,))rab{X - /x)) - (3.8) 

and this concludes our proof. □ 

The construction of the Lax pair associated to each one of the integrals of motion utilizing 
the underlying classical algebra expressed through (3.4) is illustrated in what follows. Recall 
that all the charges in involution are obtained via the transfer matrices defined in (3.3); also 
the classical equations of motion are expressed as: 

^={tH(A), L(;.)}. (3.9) 

To identify the corresponding expression for A we need to formulate the following expression: 

[tr, L:{X), L,(/.)} =ntr„ [l^-^X) {l,(A), L,(/.)}) 

= n tra{Ll-\X) [r„6(A - ^x), L„(A) + L^]) 

n tra{Ll-\X) r„ft(A - ^x)) L,{fj,) - L,{fi) n tra^L^-'iX) r„fe(A - //)). (3.10) 
Comparing (3.1) with (3.9) through (3.10) we conclude that [21, 3] 

A(A, ii)=ntra[L^^-\X) r,fe(A-/i)). (3.11) 



In the case where the r-matrix is the Yangian matrix [23] 

r(A) = \v (3.12) 

V the permutation operator V {a®h) =b® a, the expression (3.11) is modified as: 

A(A) = — — L"-i(A). (3.13) 
A — ji 

Example. Let us discuss a typical example associated to linear Poisson structures, that is 
the classical Toda chain. Consider the classical Toda chain [24] . The associated L-matrix is 
given as: 

N N-1 N-1 

L{X)=ue 2 eiN + u e 2 em + 2_^Pjejj + e 2 gj^+i + 2^ e 2 Cj+ij 

j=i j=i j=i 

(3.14) 

where we define: u — e^^, and {eij)ki — SikSji. The Lax operator (3.14) satisfies (3.4) with 
the r-matrix given by [25, 26] 

u +u 1 ^ ^ 



j=l k>] k<j 



where Ui = e 



2\i 



The momentum and Hamiltonian of the system may be now extracted. This may be 
achieved by just considering the first two generating functions trL, trL^ and expand in 
powers of u. It is convenient for our purposes here to express the L-matrix in a more 
compact form as: 

L{X) (xu^A + uB + D (3.16) 

then 

trL{X) (xu^ trA + utrB + trD (3.17) 
but due to the structure of the matrices A, B, D (3.14) we immediately conclude: 

N 

h = trB = ^Pj. (3.18) 
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To obtain the second integral of motion we expand 

trL'^iX) oc u"^ tr{A^) + tr{AB + BA) + v? tr{AD + DA + B^)+u tr{DB + BD) + tr{D'^). 

(3.19) 

Recalling the form of the involved matrices via (3.14) we conclude that the first non-trivial 
contribution, which essentially corresponds to the Hamiltonian, is given by the zero order 
terms, i.e. 

1 . N N 

72 = --tr{AD + DA + B^) = ... = -- _ ^ e9.+i-9._ (3 20) 

It is clear from the preceding discussion that by construction 

{/i, /2}=0. (3.21) 

We shall come back to the Toda model in a subsequent section providing a "dual" description 
of the model due to Sklyanin [32]-[34]. 

Exercise 3.1. Derive the exchange relations among g^, pj. 
Tip: use the following useful identity: 

^ij ^kl — ^jk ^il 

4 Quadratic Poisson structure: the discrete case 

Our main aim in this and the next section is to present a systematic means for constructing 
discrete and continuum integrable theories [1, 2]. Based on solely algebraic considerations 
we shall provide the general setting for building such systems as well as explain in detail 
how the associated charges in involution together with the corresponding Lax pairs may be 
extracted from the generating functional i.e. the transfer matrix of the system. Wc shall 
first focus on onc-dimcnsional discrete integrable models and then shall go on and deal with 
1+1 dimensional integrable field theories. Both discrete and continuum integrable theories 
share the same underlying algebra described by certain quadratic relations, which will be 
introduced below. 



(3.22) 
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4.1 Local integrals of motion 

Recall first the auxiliary linear problem in the case of discrete integrable models based on 
quadratic Poisson algebras. Lax representation of classical dynamical evolution equations [6] 
is one key ingredient in the theory of classical integrable systems [1]-[10] together with the 
associated notion of classical r- matrix [20, 21]. Introduce the Lax pair (L, A) for discrete 
integrable models [9], and the associated auxiliary problem (see e.g. [1]) 

Ipn = A„ (4.1) 

Prom the latter equations one may immediately obtain the discrete zero curvature condition: 

Ln = A„+i Ln - Ln A„. (4.2) 

The monodromy matrix arises from the first equation (4.1) (see e.g. [1, 2]) 

Ta{X) = L„^(A) Wi(A) . . . L„i(A) (4.3) 

where the index a denotes the auxiliary space, and the indices 1, . . . , denote the sites of the 
one dimensional classical discrete model. Expansion of the transfer matrix in powers of the 
spectral parameter gives rise to the charges in involution. Usually as will be transparent in 
the examples that follow the In of the transfer matrix provides the local integrals of motion. 

Consider a skew symmetric classical r-matrix, which is a solution of the classical Yang- 
Baxter equation, and let L satisfy the associated Sklyanin bracket 

{-^an(A), LbmilJ')^ = [r„b(A-/x), L„„(A) Lf,milJ')^Snm- (4-4) 

The following proposition for the monodromy matrix may be then stated. 

Proposition 4.1. Show that the monodromy matrix defined in (4.3) also satisfies the 
quadratic algebra (4.4): 

[Ta{X), n{X)} = [r„6(A-/i), T„(A) 71(/i)]. (4.5) 
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Proof. Prom the LHS of the latter equation and the definition of the monodromy matrix: 
|l„7v(A) . . .Lai(A), Lbiv(Ai) ...iv6i(At)| = 

Lan+lLbN{^) ■ ■ ■ -f'bn+l(/^)|-^an(A), Lb„(/i) |La„_i (A) . . . Lal{X)Lbn-l{^) ■ ■ ■ -^^6l(/^) = 

n 

LaN{X) ■ ■ ■ Lan+lLbNifj) ■ ■ ■ Lbn+l{l^) ^ab, Lan{X)Liya{lj) I/a„_i(A) . . . Lai (A)L{,„_i(/x) . . . I/foi (/x) — 

n 

in)... {l^)rabLan{X) . . . Lal{X)Lbn{fl) ■ ■ ■ Lfci (/x) - 

n 

LanLbwi/J') ■ ■ ■ Lbn{lj)rabLan-\{X) ■ ■ ■ Lal{X)Lbn-l{lJ') ■ ■ ■ Lfti(/x), (4.6) 

n 

but notice two consecutive terms with opposite signs cancel each other, so the only terms 
that survive in the sum are the very first and the last ones. Thus the expression above 
becomes: 

{t„(A), Tb{X)] = TabiX - ii) T„(A) Tb{ii) - T„(A) T,(/.) r„;,(A - /x), (4.7) 
which concludes our proof. □ 

Proposition 4.2. Define the transfer matrix of the discrete system as: 

t{X)=traTa{X). (4.8) 
Then show that it provides the charges in involution, that is: 

{^(A), t(A')} = 0. (4.9) 

Proof: The proof is straightforward based on the fundamental quadratic algebra defined by 
(4.5), and the definition of the transfer matrix: 

{t(A), t(A')} = [traTa{X),trbTb{>:)]=trab[Ta{X),Tb{>:)] 
= trab[rab{X-X:), T„(A)rb(A')] 

= trab[rab{X - X!) Ta{X) Tb{X') - Ta{X) Tb{X') rab{X - A')) =0(4.10) 
and this concludes our proof. □ 
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4.2 Lax pair construction 



Let us now briefly review how the Lax pair associated to each local integral of motion is 
derived via the r- matrix formulation (see also [1, 3]). Recall the Lax pair (L, A) for discrete 
integrable models, and the associated discrete auxiliary linear problem (4.1). 

First some necessary notation is introduced. We define for i > j: 

Ta{iJ;X)^Lai{X) Lai-l{X)...Laj{X). (4.11) 

Proposition 4.3. The operator A(A) of the Lax pair L, A is identified as: 

A„(A,/x) =ri(A) tr„ [T„(iV, n; A) r,b(A-/x) T,(n- 1,1; A)]. (4.12) 

Proof: To be able to construct the Lax pair we should first formulate the following Poisson 
structure [1]: 

|r„(A), Lbj{jJ,)^ = LaN{X) . . .L„„+i(A)|l„„(A), Lbn{lJ,)^Lan-l{X) . ■ ■ Lal{X) 

= LaN{X) . . . Lan+l{X) (rab{X - fl)Lan{X)Lhn{n) - Lan{X)Lbn{fl)rab{X - lx)^Lan-l{X) . . . L„i(A) 

= Ta{N, n + 1; A) r„5(A - n) T,{n, 1; A) L;,„(/x) - L;,„(/x) T{N, n; A) r„5(A - //) 7;(n - 1, 1; A). 

(4.13) 

It then immediately follows for the generating function of the local integrals of motion: 

|lnt(A), L,„(/x)} = ri(A) tr„(r„(iV,n + l;A) r„,(A-/.) T„(n,l;A)) L^(/x) 

- LM t-\X) tra(T,{N,n;X) r„(.(A - //) r„(n - 1, 1; A)) (4.14) 

Recalling the classical equation of motion 

Lnii^) = {int(A), L„(/x)} (4.15) 

and comparing with expression (4.14) we obtain expression (4.12), which concludes our proof. 
□ 

In the special case where the r-matrix is r(A) = jP [23], the operator A of the Lax pair 
becomes: 

A„(A,/x) = T(n-1,1;A) T{N,n;X). (4.16) 

A — IJ, 
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Expansion of the latter expression in powers of j (or A) gives rise to the time components 
of the Lax pairs corresponding to each local integral of motion: 



5 Examples of discrete integrable models 
5.1 The discrete NLS system or DST model 

The construction presented in the previous section will be exemplified with the study of 
the discrete version of the NLS model. This model is also know as discrete-self-trapping 
(DST) model describing the non-linear behavior of small of big molecules (see [27]-[34]). We 
introduce the associated Lax operator given by (see e.g. [30, 31]): 



-X, 1 



(5.1) 



where N ^ = 1 — x^-Xj. It also turns out, via the quadratic algebraic relation satisfied by the 
L operator, that x, X are canonical variables, i.e. 

{x„ X,} = S,,. (5.2) 

To derive the /oca/ integrals of motion one should expand the Ini(A) in powers of j. Let 
us expand the monodromy matrix: 

1 ^ 

r(A^oo) (X DN...Di + -Y^DN...Di+iAiDi_i...Di 

1=1 



i>j 

+ ^ J2 A+i A A-i ■ ■ ■ Dj+iAj . . . Dk+iAk ...Di 

i>j>k 

+ ... (5.3) 



Taking into account the latter expansion we conclude: 



In t{X ^ oo) oc jh + + + . . . (5.4) 
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where the extracted integrals of motion have the following familiar form (see also e.g. [30, 
31, 35]) 

N 
i=l 

N ^ N 

i=l i=l 

N N ^ N 

h^-Yl ^^+2X^ + J](Ni + N,+i)x,+iX, + - J] . (5.5) 

i=l i=l i=l 

The latter provide the first three integrals of motion (number of particles, momentum and 
Hamiltonian respectively) of the whole hierarchy for the NLS model. The continuum limits 
of the above quantities provide the corresponding integrals of motion of the continuum NLS 
model [30, 31, 35] as will be also transparent later in the text. It is worth noting that the 
latter expressions are valid in the quantum case as well (see e.g. [30, 31, 35]). 

The associated Lax pairs may be now identified. As already noted expansion of the 
expression (4.16) in powers of j provides the Lax pairs associated to each one of the local 
integrals of motion (see also [22]): 

^(3) ^ f /^^ + XjXj_i fXXj - XjNj + Xj+i\ 

' [-liXj^i + X,-iN,_i - X,_2 -x,X^-i J' 

Both the Lax pair via the zero curvature condition, and the Hamiltonian description 
give rise to the same equations of motion. Consider for instance the equations of motion 
associated to Is (and the Lax pair L, A^^^). Indeed from 

X,- = {73, X,}, X, = {/3, X,}, (5.7) 

and via the zero curvature condition for the pair L, we obtain the following set of 

difference equations: 

±j = Xj+2 - 2xj-+iNj - Xj+iNj+i + XjN^j + x|Xj_i + xj+i 

X,- = -X,-_2 + 2X,_iN, + X,-_iN,_i-X,-N,2_x2x^.^^_X^._^. (5.8) 

With this we conclude our brief review on the periodic discrete NLS model. 
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5.2 The Toda model 

The Toda model may be also obtained as a suitable limit of the DST model. Then the 
harmonic oscillator algebra formed by (x, X, xX) reduces to the Euchdean algebra generated 
by (e^, p) (see [32]-[34]). The associated Lax operator for the Toda model is given by: 




LajW = -7 . • (5-9) 



Exercise 5.1. Show that the first two local integrals of motion (momentum and Hamilto- 
nian) for the Toda chain are given by the following expressions: 

N 

^ N N 
j=l 3=1 

Extract also the associated equations of motion. 

Recall the linear description of the Toda chain presented in section 3; the expressions for 
the momentum and Hamiltonian clearly coincide. 



Exercise 5.2. Extract the Lax pair associated to the Hamiltonian I2 of the Toda model, 
and show that is given as 

- (J,-, t) ■ 

Show also that the corresponding equations of motion are given by: 

Pj = q^^ g^. = e*+i-* -e«^-«^-i. (5.12) 



6 Quadratic Poisson structure: the continuous case 
6.1 Local integrals of motion 

The basic notions regarding the Lax pair and the zero curvature condition for a continuous 
integrable model are reviewed following essentially [1]. Define ^ as being a solution of the 
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following set of equations (see e.g. [1]) 

— = U(a;,i,A)* (6.1) 
^ = Y{x,t,X)^ (6.2) 

U, V being in general n x n matrices with entries defined as functions of complex valued 
dynamical fields, their derivatives, and the spectral parameter A. The monodromy matrix 
from (6.1) may be written as: 

T(a;,|/,A) = Pexpj^ U(a;', t, A)(ia;'}. (6.3) 

The fact that T also satisfies equation (6.1) will be extensively used to get the relevant 
integrals of motion. Compatibility conditions of the two differential equations (6.1), (6.2) 
lead to the zero curvature condition [8]-[10] 



U- V' + 



U, V 



0, (6.4) 



giving rise to the corresponding classical equations of motion of the system under consider- 
ation. 

Hamiltonian formulation of the equations of motion is available under the r-matrix ap- 
proach. In this picture the underlying classical algebra is manifestly analogous to the quan- 
tum case. The existence of the Poisson structure for U realized by the classical r-matrix, 
satisfying the classical Yang-Baxter equation, guarantees the integrabihty of the classical 
system. Indeed, assuming that the operator U satisfies the following ultra-local form of 
Poisson brackets 

[Ua{x,X), U6(y,//)}= [r„6(A-/x), U„(x, A) + U6(y, //)] 5{x-y), (6.5) 

then T{x,y,X) satisfies (4.4), and consequently one may readily show for a system on the 
full line: 

\^lntr{T{x,y,X,)}, lntr{T{x,y, X2)}] = (6.6) 

i.e. the system is integrable, and the charges in involution -local integrals of motion- are 

obtained by expansion of the generating function lntr{T{x,y, X)}, based essentially on the 
fact that T satisfies (6.1). 
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6.2 Lax pair construction 



The construction of the time component V of the Lax pair associated to the local integrals 
of motion is outhned below. The first step is to formulate the following expression: 



|t„(L,-L,A), Ub(a:,/x)} = 
where we define 



dx 



+ 



la{x, A, jj) = Ta{L, X, A) r„6(A - jj) Ta{x, -L, A). 



(6.7) 



(6.8) 



For more details on the proof we refer the interested reader to [1], but the latter is an 
immediate consequence of the process followed in the discrete case after considering the 
continuum limit. We shall discuss in more detail the continuum limit of integrable discrete 
models in a subsequent section. 

Bearing also in mind the definition of t{X), and (6.7) it is straightforward to show: 

dY{x, A, /x) 



{in t{X), U(x,/x)} = 



dx 



+ 



Y{x,X,ii), U(x,/x) 



where we define 



Y{x,X,ii)^r\X) traMa(x,X,ii). 



(6.9) 



(6.10) 



Special examples of classical integrable field theories are presented in the subsequent 
section. 



7 Examples of integrable field theories 
7.1 The NLS model 

A particular example associated to the rational r-matrix [23] , that is the non- linear Schrodinger 
(NLS) model will be investigated below. The Lax operator for the system is given by the 
following expressions [1, 36]: 



U=M \ , (7.1) 



and it is shown via the linear algebra (6.5) that ip, ip satisfy 

{^(x), tl;{y)]=5{x-y). (7.2) 
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Our main aim is to extract the first couple of local integrals of motion via the expansion of 
In t{X) in powers of j. We shall start this construction recaUing that the classical monodromy 
matrix satisfies the fundamental quadratic algebra (4.5). Consider also the following ansatz 
for the monodromy matrix: 



T{x, y, A) = (1 + W{x, A)) e^(^'^'^) (1 + W{y, X))-\ 



(7.3) 



W and Z are purely off-diagonal and diagonal matrices respectively. We also assume that 
W, Z are expressed as: 



(7.4) 



n=l n=— 1 

Our aim henceforth is to identify the elements Wm Z^-, and hence the integrals of motion. 
It is technically convenient to split the Lax operator into a diagonal and an off-diagonal part 
as 



A /I 



V 



2 \0 -1/ VV' 



(7.5) 



Substituting the ansatz (7.3) into (6.1), and splitting the resulting equation into a diagonal 
and an off-diagonal part one obtains 

dW 



dx 
dZ^ 

dx 



+ WUd - UdW + WVaW - U„ = 0, 



(7.6) 



We may now solve the latter equations for each order and determine all W^^'^\ Z^"'^ and 
consequently the local integrals of motion. More precisely, recall the generating function of 
the local integrals of motion 

Q{X) = In (trT(A)) (7.7) 
due to the ansatz (7.3) we can substitute the monodromy matrix accordingly and obtain: 



g{X)^\ntr\{l + W{L)) e^^^'"^) {1 + W{-L))-^ 



but due to the choice of periodic boundary conditions we conclude that: 



(7.8) 



g{X) = ]ntr 



(7.9) 
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Let us first evaluate the first couple of l^^'^'s, 

\ip"{x) - \i;{x)\^i;{x) ' ■ ^ ' ' 

Similarly, through (7.6) the diagonal elements Z'-*-' are given as: 

^ ^ Z(^) = ( i^{x)i){x) 

-L) ' \ - J^^dx i){x)i}{x) ^ 

- J^j^ dx ijj'{x)i/j{x) 

— J^j^ dx il)[x)il)' [x) 




j!:^dx (^r(x)^(x)-mx)\') 



.(7.11) 



The expansion of the generating function Q in powers of j provides the local integrals 
of motion of the model under consideration. Note that due to the fact that for A — > oo the 
leading contribution comes from z['^^^ , all the local integrals of motion are given essentially 
by the quantities z[i\ More precisely the first three integrals of motion -number of particles, 
momentum and Hamiltonian - may be respectively expressed as 



Af ^ I dx ip{x)ip{x) 



1 '■^ 



dx 

L 



n ^ -J dx (I'tpix)]'' + -tp'{xy4>\x)y 



(7.12) 



The next natural step is the construction of the relevant Lax pair, in particular the asso- 
ciated V operators. Recall the generic expression of the V-operator derived in the previous 
section: 

V(x,A,//) =r^(A) tra(Tamx,X) rab{X-pi) r„(x,-L,A)). (7.13) 

Recall that in this case the relevant r-matrix is the Yangian, so substituting r(A) = y in the 
latter expression we end up: 

Y{x, A, /.) - T(x, -L, A) T(L, x, A). (7.14) 

A — fJj 
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Ultimately we wish to expand T{X) in powers of in order to determine the relevant 
Lax pairs from (7.14). Recall the ansatz for the classical monodromy matrix (7.3), after 
substituting (7.3) in the expression (7.14) and recalling that the leading contribution comes 
from the Zu term, we conclude: 

V(A) = + ^(^' ^)) ^ (1 + ^(^' ^))~' (7-15) 
A — jj, 

where we define 

-(;:)• 

We then expand expression (7.15) in powers of j, and thus we identify the corresponding 
V-operators. Indeed, 

V(A,/^) = E:X^- (7-17) 

n=l 

We provide here the quantities associated to the first three integrals of motion: 

V(i) = V, V(2)(/x) 



V(^)(/.) = r-^t ^^ + ^1, ... (7.18) 

Details on the computation of these quantities are left as an exercise to the interested reader. 

It is now straightforward to obtain the equations of motion relevant to each integral of 
motions. Focus for instance on the Hamiltonian of the system T-L (7.12). From the zero 
curvature condition for the pair U, V*^^'' one has: 



0, (7.19) 



which leads to the following familiar equations of motion 



It is clear that the same equations of motion are entailed 



kx) = ^t^-2\^l^{x)\'i^{x). (7.20) 



^(x) = [n, ^(x)}, ^{x) = [n, ^(x)}, (7.21) 
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and this of course verifies the consistency of the whole process. Results on the generahzed 
vector gl{ri) NLS model are presented in [35]. 



Exercise 7.1. Introduce the isotropic Landau-Lifshitz (classical Heisenberg) model with 
Lax operator (see e.g. [1, 37]): 

1/53 c- \ 1 

satisfying the linear algebra with r matrix being the Yangian, and thus the Poisson structure 
of the phase space for the physical quantities Si{x) is given by the Poisson brackets 

{Saix), S,{y)} = 2iSabcSc(x)5{x - y), (7.23) 

where Eabc is the totally antisymmetric Levi-Civita tensor with value £123 = 1. Based on the 
knowledge of the Lax operator U: 

(i) Determine the first two integrals of motion from the expansion of the IntrT(A) in powers 
of A: 

n=Q 

These are the momentum and Hamiltonian and are given as: 

r-L c dS2 Q dSi 



■L 1 + 



4 J _i \\dx J \ dx J \ dx 

(ii) Derive the V-operator associated to the Hamiltonian and show that is of the form: 

(iii) Show that the equations of motion relevant to the Hamiltonian are give by: 
7.2 AfRne Toda field theories 

We focus on a physically interesting class of integrable models that is the affine Toda field 
theories. The first model in this hierarchy is the well known sine-Gordon model. Suitable 
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massless limits of these theories give rise to the Liouville model and to its higher rank 
generahzations (massless ATFTs) . We shall exemphfy our study by considering the first two 
models of the hierarchy, i.e. the A^^^ (sine(sinh)-Gordon) and the A^2^ ATFT's, as well as 
the Liouville theory. 

7.2.1 The Kac-Moody ii^^ algebra 

Before we discuss the affine Toda field theories in detail it will be useful for what is described 
in the subsequent sections to recall the basic definitions regarding the Kac-Moody An^ 
algebra. Let 

Qij = 2Sij - {5i j+i + Si j_i + 5ii 5jn+i + Sin+i Sji), i, j e {1, . . . ,n + 1} (7.28) 
be the Cartan matrix of the affine Lie algebra s/n+i^ [38]. Also define: 

m 

[m]\ = 11 k, 

k=l 



m 
n 



[m\\ 
[n]\ [m - 



m > n > Q. 



n 



(7.30) 



Definition. The affine enveloping algebra sln+i has the Chevalley-Serre generators e^, /j, 
hi, i E {1, . . . ,n + 1} obeying the defining relations: 



hi 1 hj 

1 fj 

and the Serre relations 

l — Qij 

n=0 







hi 1 Gj 



Qij Cj 



hi ) fj 



5ij hi, i,j e {1, . . . ,n} 



(7.31) 



1 - aij 
n 



xl""'^"" Xj X7 = 0, e {ei: fi}: i ^ j. (7.32) 



Remark: The generators e^, fi, hi for i e {1, . . . , n} form the sln+i algebra. Also, hi — e— e^+i, 
where the elements belong to gln+i- Recall that gln+i is derived by adding to sln+i the 
elements ejie{l,...,n-|-l}so that Y17=i ^» belongs to the center. Furthermore, there exist 
the elements Sij e gln+i i 7^ j, with £i ,+1 = e^, £i+i , = /j i e {1, . . . n - 1} and 



£ij = Sik Skj - £kj ^ik, j^k^i, i, j e {1, . . . , n + 1}. 



^For the SI2 case in particular 

aij = 2Sij - 2{6ii 5j2 + ^ji), i, j e {1, 2} 



(7.33) 



(7.29) 
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It is clear that £ij E gln+i because they can be written solely in terms of the generators e^, 
fi, i e {!,..., n}. □ 

We also provide explicit expressions of the simple roots and the Cartan generators for 
An^ [39] , in order to express the Kac-Moody algebra in the Cartan- Weyl basis. The vectors 
ctj = {a} , . . . , af) are the simple roots of the Lie algebra of rank n + 1 normalized to unity 
ai • ai — 1, i.e. 
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a, - (O,... ,0,-y^,y^,0,... ,0), ^e {!,... n} (7.34) 
Also define the fundamental weights — (//^ , . . . , /jt^) , k — 1 , . . . ,n as (see, e.g., [39]). 



The extended (affine) root a„+i is provided by the relation 

n+l 



(7.35) 



(7.36) 



i=l 



The A^^ algebra in the Cartan- Weyl basis is expressed as: 

2 



ai ■ H. 



Oil 



E E 

We provide below the Cartan- Weyl generators in the defining representation: 

Eai ^ j+i , E-ai ^ ^i+l i i -^ccn+i ^ ~^n+l 1 ) -E'-a„+i ^ n+l 



(7.37) 



(7.38) 



For A2^ in particular we have: 



Q!i = (1, 0), a2 = {~,^), as 



1 

r 2 ^ 



(7.39) 



define also the following 3x3 generators 

Ei^EU^ei2, E2^Et2 = €23, Es = Et^ = -e3,, (7.40) 
where we recall the matrices Cij: {eij)ki — Sik 5ji. The diagonal Cartan generators Hi 2 are 
Hi^l (en - 622) , ^2 = ^ (en + 622 - 2633) ■ (7.41) 
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7.2.2 The A'h^ ATFT: preliminaries 



Basic notions relevant to the derivation of the local integrals of motion and the corresponding 
Lax pairs in ATFT's are reviewed. Recall first the Lax pair for a generic A^^ theory [40]-[44]: 

V(x, t,u)^^d,<^>-H+^ (u ei*-^ E+ 6"^*-^ - - e"!*-^ E_ ei*'^) 
2 4 \ u / 

iJ{x, t,u) = ^U-H+^ (u e§*-^ e-^*-« + i e'^*'^ E_ e^"") (7.42) 

u = e"+i is the multiphcative spectral parameter, and 11 are n- vector fields, with 
components (pi, TCi, i E {1, . . . ,n} which are canonical i.e. 

[Ux), 7rj{y)}^5ij5{x-y). (7.43) 

Also define: 

n+l n+1 

E+ = J2Ea„ E^ = Y,E-a, (7.44) 

i=l i=l 

ai are the simple roots, H (n- vector) and -E±a. are the algebra generators in the Cartan-Weyl 
basis corresponding to simple roots defined in (7.38). 

The associated classical r-matrix is given by [25, 26]: 

1 / •> N n+l n+l 

r(A) = Y: ^ + E e[-"(-^-)-(-)^l\., e,,. (7.45) 

smh(A) smh(A) 

Notice that the Lax pair has the following behavior: 

Y\x,t,-u-^) = Y{x,t,u), lS\x,t,u-^) =l]{x,t,u) (7.46) 

where * denotes usual transposition. 

To recover the local integrals of motion of the system under consideration we shall follow 
the quite standard procedure, and expand In trT{u) in powers of u^^. To expand the transfer 
matrix and derive the local integrals of motion we shall need the expansion of T[x, y, u). In 
what follows in the present section we basically introduce the necessary preliminaries for 
such a derivation, and later in the text we reproduce the known integrals of motion for the 
A^^\ A^^^ ATFT's on the fuU line. 

Following the logic described in [1] for the sine-Gordon model, we aim at expressing 
the part associated to £'+ in U, independently of the fields, after applying a suitable gauge 
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transformation. More precisely, consider the following gauge transformation: 

T{x, y, u) = n{x) f{x, y, u) fl'^y), n{x) = e^^^^^'^. (7.47) 

Then from equation (6.1) we obtain that the gauge transformed operator U may be expressed 
as: 

u(x, t, u) = n-^ {x) u{x, t, u) n{x) - n-^ {x) (7.43) 

After implementing the gauge transformations the operator U take the following simple form: 

U(x, t,u)^^e-H+j (uE+ + , (7.49) 

where we define: 

e = n - x_ = e-^^-" e_ e^^-", (7.50) 

© is a n- vector with components 9i, and it is clear that T now satisfies: 

^^^^^^^ = U(x,A)f(x,y,A). (7.51) 

Consider now the following ansatz for T as — >■ oo [1] 

f{x, y, u)^{I + W{x, u)) exp[Z(x, y, u)] (I + W{y, u))-\ (7.52) 

where W, W are off diagonal matrices i.e. W = ^j^^ WijEij, and Z, Z are purely diagonal 
^ = S"=i^ ZiiEii. Also 

fc=-l k=0 

Inserting the latter expressions (7.53) in (7.51) one may identify the coefficients wj^^^ and 
Z^^\ Indeed from (6.1) we obtain the following fundamental relations: 

— = u(^) + (u(^) PF)(^) 

dx 

^ + Wi](-D) - U(^)iy + iy(u(«)iy)(^) - _ (u(0)^)(0) = 0, (7.54) 
dx 

where the superscripts O, D denote off-diagonal and diagonal part respectively. We shall 
come back to the solution of the latter set of equations in the subsequent sections when 
investigating particular physical examples. 
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7.2.3 The sine-Gordon and Liouville models 



Note that we focus here on the sine-Gordon and not the sinh-Gordon model, as one would 
expect from the notation of the previous section, where we introduce real ATFT's. In 
fact the only difference compared to the notation introduced previously is that (3 ^ if3; 
this is basically the only modification occurring between real and imaginary ATFT's. The 
associated classical r-matrix [25, 26] in this case is expressed as 

1 /(^) cosh A (J- \ . . 

and the Lax operator for the sine- Gordon is now given as 

V{x,t,u) = ^n{x)a^ + I^e^^-V^e"^^'^^ - ^e-^^^V^e^^-^ (7.56) 
^ 4i ^ At At ^ ^ 

g.x,y,z ^Yie familiar 2-dimensional Pauli matrices, u = e^. U satisfies the fundamental 

quadratic algebraic relation with the classical r-matrix (7.55), which gives rise to: 

[<l>{x), 7T{y)}^5{x-y). (7.57) 

To consider the formal series expansion of T we shall need the following symmetry of the 
Lax operator: 



and consequently 



U(ii-\0,7r) =U(-u. -0,7r) (7.58) 



r(xi"\0,7r) =r(-w,-0,7r). (7.59) 



As in the previous section we aim at expressing the term of order u in V independently of 
the fields, after applying a suitable gauge transformation [1]. More precisely, consider the 
gauge transformation (7.47) with 

= e'^^'^^ (7.60) 

then the gauge transformed operator U is expressed as: 

U(x, t, u) = ^f(x)a^ + - I^e-^-^-V^ef '^'^^ (7.61) 

4z Ai Ai 

where we define 



f(a;,t) = Ti{x,t) + (p'{x,t). 
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(7.62) 



We consider the decomposition (7.52) for T, as — > oo [1]. W, is off diagonal matrix 
and Z is purely diagonal, and are expressed as in (7.53). Inserting expressions (7.53) in 
(7.51) one may identify the matrices W^''^ and Z^'^\ via (7.54). 

It is sufficient for our purposes here to identify only the first couple of terms of the 
expansions. Indeed based on equation (7.54) we conclude (see also [1]): 



m 

y^(2) ^'^^y _ ^sin(/30) (jy-^a^ 



m 



2im? 



(7.63) 



Then, as in the previous example we also identify the diagonal elements In particular 
from the first of the equations (7.54) we extract the following expressions {x — L, y — —L): 



2 



e 



(7.64) 

Notice again that for —iu^oo the leading contribution comes from the Zu elements, due 
to the form of Z^~^\ see also similar argument for the NLS model described previously. 

Then recalling the expression for the generating function of integrals of motion, intro- 
duced in the previous example, we conclude that 

,(1) ^ ^(1) ^ £ ( - + cos^0) . (7.65) 

If we now do the same kind of expansion but for A — > — oo we basically end up to a similar 
expression as before, by simply exploiting the basic symmetry of the monodromy matrix 



T{u-\(f), 7r) = T(-u, -(f), tt) (7.66) 



we find that 



where we define 

f{<f>, 7r) = f(-0, tt). (7.68) 
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In fact it is clear that any combination of the quantities is also one of the charges 

in involution. So we find that the Hamiltonian is 

I^'^ + I^-'^o,U^ j\x (i(7r2 + 0'2)_ ^cos/30), (7.69) 

and we may also identify the momentum as: 

/(I) _ /(-I) oiV^ J dx (t)'{x) 7r(x). (7.70) 

Let us also identify the V operator associated to the Hamiltonian and momentum of the 
sine-Gordon model. Recall the generic expression for V 

V(x,A,/x,x) =r^(A)tr„(r„(L,x,A) r„6(A - /x) r„(x,-L,A)). (7.71) 

Recall the ansatz for T (7.52), and that the leading contribution comes from the Zn, then 
the latter formula may be expressed as: 

V(x, A, ix) = [(1 + Wa{x))-^ n-\x) r„,(A - //) n{x) (1 + W^ix))] (7.72) 

Our aim now is to expand the latter expression in powers of (recall u — e^). The 
r-matrix is expanded as 

fl-\x) r(A - fi) n{x) = f(°) + + 0{u-^) (7.73) 

and 

1 + W{x) = A + u-^W^^^ + 0{u~^) 
(1 + W{x))-^ = A-^ + u-^f^^ +0{u-^), (7.74) 



where we define: 



2 \^ -a' + lj \e2(T+ 



V — e^. 



We keep for our purposes here only the first terms of the expansion; these arc sufficient to 
provide the relevant V quantities. Gathering all the information from the latter expression 
we eventually obtain: 

V(i)(x,/i) = ^f{x,t)a' + vn{x,t)ayn-\x,t). (7.76) 

^ III 
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Similar computation for the expansion in powers of u can be made and the final result in his 
case is directly given as: 

Y^-^\x, n) = -7^f(x, t)a' + v-^n-\x, t)ayn(x, t). (7.77) 

By adding the two expressions above we end to the desired result (we have also multiplied 
with an overall factor ^) 

V(«) = t)a' + "^^{x, t)am-\x, t) + t)ayn{x, t). (7.78) 

Similarly, subtraction of the two quantities Y^^\ V(^^) leads to the time component of the 
Lax pair associated to the momentum, which coincides with U. 

Note that appropriate massless limit of the sine-Gordon model gives rise to the Liouville 
theory through a particular limiting process (see e.g. [45, 46]). The entailed Lax pair of the 
Liouville theory takes the form: 

^^^^ if -in{x) \ ^^^^ I ( 2e-^-*'^(^) 

2 y 4sinh(A — i0(x)) i7r(x) y ' 2 I 4cosh(A — i0(a;)) i(f)'{x) 

The Lax pair satisfies the zero curvature condition, which leads to the corresponding equa- 
tions of motion i.e. 

sine-Gordon model: (f){x) — 4)"{x) -\- — sin(;50(x)) = 

Liouville model: 0(x) - 0"(x) - 4ie-2^^(^) = 0. (7.80) 



Exercise 7.2. Starting from the U operator of the Liouville model derive the energy and mo- 
mentum of the model well as the associated V operators. Determine also the corresponding 
equations of motion. 

7.2.4 The A^^^ ATFT 

We shall now study the second member of the ATFT hierarchy, i.e. the model. It will 
be useful in what follows to introduce some notation (see also [49, 50]): 

|e • = diag(a, 6, c), f © • ^ = diag(a, 6, c), e'^"-* = 7^ (7.81) 
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©($, n) = ©(— n) explicit expression of a, 6, c and 7j can be found in Appendix A. 
Prom the first of equations (7.54) we may derive the matrices Z. Indeed one may easily 
show that: 

v^3? + C6 =0 

V -w^i? + cc y 



dx 



m 

T 



T 



v 



0. 



(7.82) 



W^2? + CcJ 



It is clear that the latter quantities are zero because of the form of wj^^ , wj^^ see Appendix 
A. Also the higher order Z^''\ Z^^^ are given by: 



dZ^ 
dx 

dZ^ 
dx 

k>0. 



m 

T 



m 

T 



73 I 
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(7.83) 



The computation of W, W is essential for defining the diagonal elements. First it is important 
to discuss the leading contribution of the above quantities as \u\ oo. To achieve this we 
shall need the explicit form of Z^~^\ Z^~^^: 



Z^-'\x,y) 



m{x — y) 



( 



e 3 



e 3 



m{x — y) 



Z7I 

e 3 



e 3 



■1/ 



(7.84) 



The information above will be extensively used subsequently. Prom the formulas (7.54), 
(7.83) the matrices W^^\ Z^^\ Z^^^ may be determined (see Appendix for details). 

The local local integrals of motion in the periodic case, emerge as usually from the 
expansion {\u\ — >■ oo) of: 

In [trT{u)] = In [ir|(l + W{L,u)) e^^^'-^'") (1 + W{-L,u))-^^ . (7.85) 
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Notice that in the case of periodic boundary conditions we put our system in the 'whole' 
hne {x — L, y — —L), and consider Schwartz boundary conditions, i.e. the fields and their 
derivatives vanish at the end points ±L. Bearing in mind that as it — > — oo the leading 
contribution of e^, (e^) (see (7.84)) comes from the e^^^, (e^^^) term, the expression above 
becomes 

In [trT(« ^ -oo)] = ^ (7.86) 

k=-l 

To reproduce the familiar local integrals of motion we shall need both Z{L, —L, u), Z{L, —L, u) 
Let 

P J-L P -^-^ 

I^-'^--'^zi^{L,-L,u)^ r^-(E^l + ^Ee^-^) (7.87) 

the momentum and Hamiltonian (and the higher conserved quantities) of the ATFT are 
given by: 

V, = J(/(-^) - /(^)) = Tdxf^ <t>'r - < <t>) 

^ i=l ^ ^ 

= ^(/« + /(-)) = /' dx(y^{.l + 0?) + ^ E ■ 
^ -^-I^ ^ i=l p i=l ' 

(7.88) 

Note that the boundary terms are absent in the expressions above, since we considered 
Schwartz type boundary conditions at ±L. Also, in the generic situation, for any ^ the 
sum in the momentum V\ and the kinetic term of the Hamiltonian "Hi runs from 1 to n, 
whereas the sum in the potential term of the Hamiltonian runs from 1 to n + 1. A study of 
generic integrable boundary conditions in ATFT's is presented in [47]-[50]. 

Exercise 7.2. Based on the process described for the derivation of Lax pairs, determine 
the V-operator associated to the integral of motion T^i (7.88). 
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8 The continuum integrable limit 



The main aim in this section is to introduce a consistent continuum hmit of integrable 
discrete models in such a way that integrability is preserved (see also [51]). The starting 
point is the lattice monodromy matrix, which is expressed as: 

Ta — LaN LaN-l ■ ■ ■ Lai- (8.1) 

Assume that L admits an expansion in powers of S as 

Lai = l + SVai + C»(5') , (8.2) 

then consider the product: 

N oo 

i; = n(l + 5U„, + 5^5"Ul?)). (8.3) 

i=l n=2 

Expanding the expression above in powers of 6, we get 

i i<j i 

These, multiple in general, infinite series of the products of local terms, are characterized by 
two indices: the overall power n of 6, and the number m of the set of indices i over which 
the scries is summed. Note that, in the T expansion one always has m. The continuum 
limit soon to be defined more precisely, will entail the limit 5 — )■ with 0{N) = 0{l/6). We 
now formulate the following power- counting rule, that is terms of the form (see also [51]) 

m 

J2 J2^,-n, (8.5) 

i\<i2<---im. 3=1 

with n> m are omitted in the continuum limit. The latter is defined by 

5Y,^ai^ I dx\}a{x) (8.6) 

■ J-L 

and similarly for multiple integrals. Here 2L is the length of the continuous interval defined 
as the limit of N5. In other words, contributions to the continuum limit may only come 
from the terms with n = m for which the power 5" can be exactly matched by the "scale" 
factor N'^ of the m-multiple sum over m indices i. In particular, only terms of order one in 
the 5 expansion of local classical matrices Lai will contribute to the continuum limit. Any 
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other contribution acquires a scale factor S'"~"^ — > 0, when the continuum hmit is taken (see 
also [51]). This argument is of course valid term by term in the double expansion, and it 
always has to be checked for consistency. 

The continuous limit of T, hereafter denoted T, is then immediately identified as the 
path-ordered exponential from y — —L to x — L 

r^Pexp([ dxV{x)] , (8.7) 



where suitable (quasi) periodicity conditions on the continuous variables of the classical 
matrix U(x) are assumed. Of course the derivation of a continuous limit requires that the 
L-matrices are not too inhomogeneous (e.g. L-matrices at neighbor sites should not be too 
different. 

The above identification of T has been built so that to straightforwardly generate the 
classical continuous limit of the Hamiltonians from the analytic expansion 

oo 

tr(r(A)) = J](A - Ao)"?^(") . (8.8) 

n=l 

We thus characterize U{x) as a local Lax matrix yielding the hierarchy of continuous Hamil- 
tonians 'H^^\ In order for this statement to agree with the key assumption of preservation 
of integrability we are now lead to require a Poisson structure for U compatible with the 
demand of classical integrability of the continuous Hamiltonians. Indeed, such a structure 
is deduced as the ultra-local Poisson bracket 

|Ui(x,Ai), U2(?/,A2)}= [ri2(Ai-A2), Vi{x, X,) + U2{y, X2)\s{x - y) , (8.9) 

where r is the classical matrix characterizing the exchange algebra of the L-operators. More 
specifically, recalling that Lai = 1 + SI}ai + 0{6'^), plugging it into (4.4) and assuming ultra- 
locality of Poisson brackets one gets 

{u„,(Ai), U5,-(A2)} = [rabiXi-X^), U„,(Ai)+U;.,(A2)]^ . (8.10) 

One then identifies, in the continuum limit 5 — )■ 0, the factor Sij/S with 6{x—y). Reciprocally, 
it is a well known result (see, for instance [1]) that if U(x) has a such an ultra-local linear 
Poisson bracket (8.9) the full monodromy matrix between —L and L has the quadratic 
Poisson bracket structure, thereby guaranteeing Poisson commutation of the Hamiltonians. 

We now focus on the continuum limit of the zero curvature condition, and recover the 
continuum expression as a suitable continuum limit of the discrete one. Moreover, the generic 
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continuum expression for the V-operator is also ensued via this process. Recall the discrete 
zero curvature condition: 

L„(A) = A„+i(A)L„(A) - L„(A)A„(A). (8.11) 
Consider the following identifications: 

Ln^l + 5U(x), An^Y(x), An+i^Y(x + 5) (8.12) 

then expressing 

Y{x + 6)= Y{x) + 6V{x) + 0{6'^) (8.13) 
and keeping the first non trivial contribution which is of order 5 we conclude: 

U(A, x) - V'(A, x) + [U(A, x), V(A, x)] = 0, (8.14) 
which is nothing else but the familiar continuum zero curvature condition. 

Example. Let us work out a particular example to demonstrate how this process may 
be applied. Consider the discrete NLS model; first introduce the spacing parameter S in the 
Lax operator (see also [52] and references therein): 

we have essentially rescaled the fields x, X and the spectral parameter as: 

X^SX, x^5x, X^5X. (8.16) 
After the aforementioned rescaling the first local integrals of motion become: 

3 

N-l N 

I, = -S' J2 - 2 E(l - . (8.17) 

j=i j=i 

We shall consider the continuum limit of the first two integrals just to illustrate how this 
works. Bear in mind the following identifications: 

Xj — >■ x(x), Xj — >■ X{x) 

Xj+i x{x + 6), Xj+i X{x + 5), (8.18) 
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also consider the approximation for the sum 

It is then straightforward to show that the first integral of motion, after keeping the first 
non trivial contributions of order 5, becomes in the continuum limit: 



J x{x)X{x). (8.20) 



Similarly, the second integral becomes after taking the continuum limit and keeping the next 
nontrivial order 5^, 

X2 = ^/" dx (x'{x)X{x) -x{x)X'{x)y (8.21) 

J — L 

It is also clear that the continuum limit of the associated Lax operator (8.15) takes the 
desired form: 

L{X)^l + 6l]{x) + 0{S^), (8.22) 

where we define 

notice that U may be always derived up to an additive constant. Take also the continuum 
limit of the A operator associated to the second charge H2 this then becomes via (8.16), 
(8.18) 

A„^VW, VW ^ J^^j ='1,^'] . (8.24) 

By identifying 

x{x) = ipix), X{x) = -il){x) (8.25) 

we end up to the familiar continuum expressions extracted in the previous section. In par- 
ticular, Ii and I2 are proportional to Af and V respectively. This is a strong consistency 
check verifying the validity of the applied continuum limit process. 

Exercise 8.1. Determine the continuum limit of the discrete NLS Hamiltonian I-?. 
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A W and Z matrices for the ^ ATFT 



We write below explicit expressions of the Z*^") matrices for the first orders: 



m 



v 





/ 


e- 


1 










e 3 


-1 










2in in 











I 


6 3 6 3 


J 









2in \ 










e 3 a c ' 
















m 


—a 


b 








esc 


-b j 






4 




(A.l) 



The higher order quantities are more complicated and we give the matrix entries below for 



and 



(2) 
12 



21 



(2) 
13 



w. 



(2) 
31 



(2) 



23 



'32 



w 



(2) 
13 



^(-273 + 71 + 72) + y (2a' + 6') + ^(-20^ - 6c), 



6 3 



-273 + 71 + 72) + 



— - — (a' - d) + (c^ - ah) 



^(-272 + 71 + 73) + y (-&' + c) + ^(6' - ac), 
1(2^2 - 71 - 73) + y (-«' - 2c') + ^(2c2 + a6), 
-^(271 - 72 - 73) + ^(26' + c') + ^(-26^ - ac) 



6 3 (^e" 
■ — (271-72-73) + — ^ 



2 

-a' + 6') + ^^(a2-6c) (A.2) 



6 3 



;-272 + 7i + 73) + 



3 3 



(c2-a6), 



^(-272 + 71 + 73) + j(26' + a') + ^(-26^ - ac) 
-^(-271 + 73 + 72) - y (2a' + c') + ^(2a2 + bc\ 



e 3 



(271 - 72 - 73) 



a2 — 



(6'-a') + ^ {-V + ac) 



= -^(-273 + 72 + 71) + y (a' - c') + ^(-a^ + 6c) 
3? = ^(-273 + 71 + 72) + y (^' + 2c ) + y (-2c^ - ab) 



(A.3) 



35 



where the prime denotes derivative with respect to x, also a, b, c, and 7^ are defined in 
(7.81) and have the following exphcit forms: 

Moreover using the expressions above and (7.83) we have: 



3 



a/j\{ earn, . ^ / , Ce - , 72 2\ 

7^(1) ilL /- ijL /- ijL 

^ = ^^(7. + 7. + 73) + V(^' - «') + + ^ + 



3 



= -t(^' + + - — (" - " ) + — (" + " + 

^ = -~(7i + 72 + 73) + |(a' - e') - ^(a^ + » + (A.5) 
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